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Abstract —We construct two-dimensional codes for correcting 
burst errors using the finite field Fourier transform. The encoding 
procedure is performed in the transformed domain using the 
conjugacy property of the finite field Fourier transform. The 
decoding procedure is also done in the transformed domain. Our 
code is capable of correcting multiple non-overlapping occurrence 
of different error patterns from a finite set of predefined error 
patterns. The code construction is useful for encoding data in 
two dimensions for application in data storage and bar codes. 

Index Terms —2D error correcting codes, 2D finite field Fourier 
transform, cyclic codes. 

I. Introduction 

Data storage technologies in magnetic and flash memories 
are shifting towards a two-dimensional (2D) paradigm since it 
offers better signal-to-noise (SNR) ratio and format efficiency 
m, inturn leading to higher storage densities. Examples of 
such storage technologies include 2D magnetic recording 13 , 
bit-patterned media 0, flash storage flu, optical holographic 
recording 0 etc. The native format of data applicable to these 
technologies is inherently two-dimensional. This is useful 
for combating 2D intersymbol-interference, noise and other 
impairments. Storage channels have noise bursts along with 
an intersperse of random errors 0. The data to be encoded 
and written onto the storage medium must be capable of over¬ 
coming a mixture of burst and random errors. Thus, designing 
efficient 2D error correcting codes along with modulation 
codes is important for these storage channels along with so¬ 
phisticated signal processing for timing recovery, equalization 
and detection 0 as a precursor step before decoding. In this 
paper, we develop a theory for encoding 2D data resilient to 
burst errors. 

While there are many well known approaches for correcting 
ID burst errors, the design of codes for correcting 2D errors 
is non trivial for the following reasons: (a) The shape of 
a 2D burst can be arbitrary, unlike a ID burst string of 
errors, (b) The errors can be correlated in 2D and non- 
separable, requiring a 2D frame work for efficient encoding 
and decoding. In the following sub-sections, we will survey 
2D coding techniques relevant within the context of our work 
and highlight our contribution. 

A. Prior work and context 

The fundamental theory of 2D cyclic codes was first for¬ 
mulated by Imai 0 with the construction of generator and 


parity check tensors for 2D codewords. The construction was 
analogous to the construction of ID linear block codes. This 
code could correct only a single burst error of an arbitrary 
error pattern and was not designed for correcting single burst 
of more than one error pattern. In a recent paper by Yoon and 
Moon 0, they provided the construction of a code with dis¬ 
joint syndrome sets among specific error patterns, improving 
the construction in 0. This helped in correction of a single 
occurrence of an error pattern from a finite set of different 
error patterns. However, they did not address the correction 
of multiple burst errors. Others Ea-DH have constructed 
codes for correcting a small cluster of errors by converting the 
2D binary codewords into symbols of non-binary codewords 
by grouping the bits along rows and columns. Although 
they can correct an error pattern of any shape, their method 
could not correct multiple disjoint cluster of errors. Blahut 
lfl2l introduced higher dimensional finite field transforms for 
realizing algebraic codes over curves. Madhusudhana and 
Siddiqui 03 extended Blahut’s decoding algorithm for 2D 
BCH codes. Analogous to the ID BCH code, the 2D extension 
is required to have a block of contiguous positions set to 
zero in the frequency domain, where, the syndromes are 
calculated sequentially. They show correction of random and 
burst errors from a deconvolution viewpoint. However, their 
work does not incorporate the idea of disjoint syndrome sets 
which otherwise can lead to miscorrections. Shiozaki DU 
presented a new class of 2D codes over GF(q) using 2D 
Fourier transform techniques of dimension (q — l) 2 . In that 
work, an exact relationship between the parity check symbols 
and the minimum distance was also provided. However, the 
code construction does not accommodate disjoint syndrome 
sets. Also, the details of the decoding procedure are not 
explicitly provided. In this paper, we present a construction of 
2D codewords in the frequency domain for correcting multiple 
burst errors. Our work is based on the 2D finite field transform 
inspired from the work in lfl2ll . 

B. Our approach and contribution 

We present a construction of 2D linear codes using the 
finite field Fourier transform for correcting multiple non¬ 
overlapping occurrences of different error patterns from a 
finite set of predefined error patterns. Our encoding scheme 
uses the conjugacy property of the 2D finite field Fourier 
transform. Unlike prior approaches 0, the encoding process 
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is simple as it does not need any parity check or generator 
tensors. Further, our code construction is tailored for correcting 
a set of predefined error patterns as well as an arbitrary 
linear combination from the set of predefined error patterns. 
Most ID decoding approaches follow a two step process of 
identifying error locations followed by error correction. In the 
2D case, we first need to identify which error pattern has 
occurred. With this information, we try to solve the erroneous 
locations. In case we do not identify the error pattern, we 
declare an occurrence of uncorrectable error. In our work, we 
also consider the correction of multiple disjoint bursts of the 
same pattern assuming the number of bits to be corrected falls 
within the minimum distance bound, a feature not considered 
in m-m. The set of common roots among all the codewords 
in the code forms the set of common zeros 0. To facilitate 
the decoding procedure, we formulate the common zero set 
accordingly. A subset of the common zero set is defined as 
the set of indicator zeros. This set of zeros are used to identify 
the error pattern that might have occurred. The remaining 
set of zeros are used to form the syndrome equations in the 
transformed domain for solving the erroneous locations in 
contrast to the approach in 03 , 03 - The entire common 
zero set is also chosen in a way to obtain disjoint syndrome 
sets among a set of predefined error patterns/error events. Yoon 
and Moon 0 gave a criterion to obtain disjoint syndrome sets 
among specific error patterns. We interpret the same criterion 
from a transform domain perspective. The construction in lfj~3t 
required a block of contiguous zeros in the transform domain. 
However, in our work, we relax this constraint and choose the 
zeros for a-priori identification of error patterns, and obtaining 
disjoint syndrome sets among predefined error patterns. 

The paper is organized as follows. In Section II, we 
highlight the code description and motivate the 2D finite 
field Fourier transform (FFFT). In Section III, we develop 
an encoding scheme for mapping 2D binary arrays into 2D 
codewords using FFFT. In Section IV, we discuss the decoding 
of 2D codewords in the transform domain that can correct 
simultaneous occurrence of disjoint error patterns with illus¬ 
trative examples. We conclude the paper in Section V. 

II. TWO DIMENSIONAL BINARY CYCLIC CODES 


where, a{x,y) and b(x,y) are polynomials over GF{ 2). 
Consider a bi-variate polynomial f(x,y). A point (x',y r ) is 
said to be a root of f(x,y) if f(x',y') = 0. For a code of 
dimension n x to, let 7 and (3 be the primitive n th and m th 
roots of unity respectively i.e., 7" = 1 and 3 rn = 1. All roots 
constitute the set V given by V = { (7*, /3 J ) |0 < * < n—1, 0 < 
j < m — 1}. The roots common to all the codewords in the 
code are referred to as the set of common zeros. For any set of 
zeros A, common to some arbitrary polynomials over GF{ 2), 
if (S, a) € A then ^i5 2 ^, er 2 ^ £ A for k = 1,2,..., n — 1 and 

n is the least positive integer for which (<5, cr) = ((5 2 , cr 2 ). 
These points are known as conjugate points in A. From the 
set V, we select a subset V c , which will be our set of common 
zeros. A 2D code is completely characterized by this set of 
common zeros. 


B. 2D Finite Field Fourier Transform (FFFT) 

The 2D finite field Fourier transform (FFFT) for an array 
of dimension n x to having elements from GF(2) is defined 
as 

n—1 m— 1 

Cg,* = EE l i6 P j * c i,r ( 3 ) 

2=0 j = 0 

Cg^ £ GF(2 X ), 2 X = lcm(n,?n) + 1, where lcm(.) 
denotes the least common multiple. If we have n = 3 and 
to = 5, the elements of Cg $ are from GF{ 2 4 ). Thus, a 
binary codeword gets converted to a non-binary codeword of 
the same dimension. This is a direct extension from the ID 
FFFT defined in 02 


C. 2D Finite Field Inverse Fourier transform 

The 2D finite field inverse finite field Fourier transform 
(IFFFT) is defined as 




1 1 

mod(n,p) mod (to, p) 


n—1 m— 1 


E E r^p-^cg^ 

0=0 0=0 


(4) 


where, p is the characteristic of the field and mod(n,p) is 
defined as the remainder obtained while n\p. 


A. Code Description 

A 2D binary code of dimension n x to is a collection 
of arrays having n rows and to columns with elements from 
GF(2). A 2D array can be described using a bi-variate 
polynomial expression as, 

c{x,y)= E °i,jX l y J , (1) 

(ij)eo 

where, f l = {(*, j)|0 < i < n — 1,0 < j < m — 1} and 
c itj £ GF( 2). 

A 2D code is cyclic if x l y : ’c(x, y) £ ‘if, Vc(ir, y) £ c &. A 2D 
binary code is said to be linear if Ci(x, y) £ ‘if, Vci(x, y) £ 
Y . Let )f [H] denote the set of all bi-variate polynomials. If 
f(x,y) € then 

{f(x,y)}n = f(x,y)+a(x,y){x n -l)+b(x,y)(y m -l), (2) 


III. Encoding 

We highlight the encoding procedure in 0 in this section 
for completeness. The encoding of the 2D codewords in [8| 
begins by selecting the common zero set V c which is followed 
by the creation of the parity check tensor. The generator tensor 
is then created from the parity check tensor and used for 
encoding a message. Later in the section, we describe the 
encoding process of the code using 2D FFFT, and highlight 
its various advantages over 0. 

A. Parity Check Tensor 

The parity check tensor 0 can have two forms. In the first 
form, each element of the parity check tensor is a binary vector 
of length equal to the cardinality of the common zero set. The 
other form has a bi-variate polynomial as its elements. The 
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second form is required because we need to generate our 2D 
codewords in polynomial form. 

Two zeros of V c are said to be equivalent if the first 
component of their coordinate are conjugates over GF{ 2). 
Thus, to find the parity position, we have to partition V c 
into equivalence classes. Two roots (e, rf) and ( 0 , y) are said 
to be equivalent © i-e-, (e,r?) ~ (0,/z), for some positive 
integer k, if e 2 = 9. Each component of the parity check 
tensor, which is a binary vector, is split in length depending 
on the size of each of the equivalence classes. For a particular 
(i,j) coordinate of the tensor, we express 7 J /3 J in terms of 
GF (2l Vc fcl) for the k th segment. The binary vectors created 
at the parity positions are all linearly independent. The binary 
vectors at the message positions can be expressed as a linear 
combination of these vectors. Let the set of all coordinates 
of the 2D code be defined as f2, and, the set of all parity 
positions be defined by the set 7r. Each vector at position (i,j) 
is expressed as 

hi,j = E 

(k,l)£n 

where, h^’P £ GF( 2) and (i. j) £ fi. To get the algebraic 
form of the parity check tensor, we express the components 
in bi-variate polynomial form as 

hij(x,y)= ^2 hppx k y l . (6) 

(/c,Z)G 7T 

The operators T x and T y (8) operating on h, :) (x, y) is 
defined as 

Tx[hi,j i. x i t/)] ^mod((i+l),n),j(^) 2/)) (T) 

■fy [hi,j ( x i 2/)] = ^i,mod((j-|-l),m) 1/) ■ (&) 

From the construction of the parity check tensor in the 
algebraic form, we see that, ho,o( x ,y) component is always 
1. Few properties of T x and T y a are as follows: 

1) T*T 3 = h itj (x,y). 

2) f(T x ,T y )= J2 fijhijix.y). 

A residue of a polynomial f(x,y) is defined as f(T x T y ). A 
bi-variate codeword f{x,y ) is said to be valid if its residue is 
zero, i.e., 

f(T x ,T y ) = 0. (9) 

B. Generator Tensor 

The relation between ( i,j) components of the generator 
tensor and the parity check tensor is given by, 

9i,j( x ^y) = ( 10 ) 

From equation ([9]), it is easy to verify that each component of 
the generator tensor is a valid 2D codeword. 

Let the message polynomial be 

m(x, y) = ^2 m i,jX l y 3 , 
where rriij = 0 for (i,j) £ it. 


The work done in a does not provide an explicit math¬ 
ematical expression for encoding. We prove a result on the 
encoding of a message vector m(x,y) into a valid codeword 
c(x, y) in the following Lemma. 

Lemma 1: The encoding equation of the message vector 
m(x, y) into a codeword c(x, y) follows as, 

c(x,y)= ^2 m i,i9i,A x iV)- (H) 

Proof: We know that a codeword c(x. y) is valid if 
c{T x ,T y ) = 0. So, we can write the codeword c{x,y) as 

c(x,y) = m(x,y) — m[T x T y ). (12) 

Simplifying further by using properties of T x and T y , and 
using (TO} in (12} , we have 

c(x,y) = J2 m iJ xl y 3 - m Li T x T P 
= m i,j9i,j(x,y), 

thereby, proving the result. ■ 

C. Finite field transform of a code 

The original form and the transformed representation of 
codewords are defined in the time and frequency domains 
respectively. We use the FFFT technique to transform the 
codewords in the frequency domain. The data is encoded 
directly into the transformed domain. This procedure is easier 
because it helps us to save space for not storing the parity 
check and generator tensors, which are of considerable size. 
This reduces the space complexity at the transmitter side 
extensively. The binary codewords are obtained from the trans¬ 
formed domain using 2D IFFFT. In the transformed domain, 
the conjugacy constraints are used to express the coordinates 
of the transformed codeword. 

Lemma 2: For a 2D code over GF(q), 

Qmod(^,n}),(mod(q0,m))' (13) 

Proof: For an element k £ GF{q), we have k q = 
k. Also, we know in case of finite fields, for elements 
{cii, 02 , ■ • • , a„} £ GF(q), the following relation holds: 

(Oi + 02 + ’ ’ ’ + <ln) q = ( a i + O 2 + ’ ’ ’ + off) . 

Thus, from (3}, 

N-l M — l 

C b = E E i iq 6 p jq ^lp 

i =0 j =0 

N-l M—l 

= E E 7 i( «v w) cij, 

i =0 j =0 

C (mod(q6 ,n)) ,(m.od(q(f),m )). 

The last step follows since 

mod (qO , n) = qO , 
mod(g<^, to) = q<j>, (14) 

thereby, proving the result. ■ 
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D. Validity of a codeword in the frequency domain 

Lemma 3: Let ( 7 9 ,/3^) be a zero of the codespace. In the 
frequency domain, Cg^ = 0 . 

Proof: We have, 

c(x,y)= E 

(i,j)e fi 

Plugging in (ff 9 ,fi^) for (x,y), we have 

E W*?* = c(y e ,f3*), 

or, E daft* ft* = °- (15) 

From Q, Cg^ = 0, which proves our result. ■ 

Thus, the condition for checking the validity of a codeword 
in the frequency domain is 

C e ,4, = 0 , (16) 

f(0,4>), where G V c . 

Example 1: Let the common zero set in the frequency 
domain be 

V^ c = {(0,0), (1,1), (1,4), (2,2), (2,3)}. 

From Lemma 3, the components in VT r will be zero. We 
show one valid codeword in the transform domain that satisfies 
Lemma 3 and the conjugacy constraints. As an example, 

/ 0 a 13 a 11 a 14 a 7 \ 

C = a 10 0 a 6 a 9 0 . 

\ a 5 a 3 0 0 a 12 ) 

Now, using the 2D IFFFT equation from ([4]), we have the 
time domain codeword as 

/ 1 1 1 0 0 \ 

c = I 0 1 0 0 0 . 

\ 0 0 0 0 0 / 

Our encoding scheme clearly avoids the creation of parity 
check and generator tensors to encode a message. Given the 
considerable size of these matrices it helps to save space for 
storing them. Lemmas 2 and 3 are sufficient to encode the 
code in the frequency domain. However, codewords are always 
stored in binary format in storage devices. This is where the 
2D IFFFT operation comes into play. After encoding, the 
codeword undergoes the inverse transformation and stored as 
binary data. 

IV. Decoding 

In this section, we will discuss decoding in the frequency 
domain. In time domain, the received codeword r, transmitted 
codeword c, and the error vector e are related by the following 
equation: 

~ “I" ) (17) 

where, 0 < i < n — 1 and 0 < j < m — 1. As FFFT 
is an linear operation, the same equality holds. Let R be the 


received codeword, C be the transmitted codeword, and E be 
the error vector in frequency domain. We have. 


{7?}e,0 — {C}e ,0 + {E}g^. 

(18) 

According to Lemma 3, 


{R\e,<t> = {E}g^, 

(19) 


V ($,(/>) where, (y 9 ,^) G V c . 

When we receive a codeword, the coordinates corresponding 
to the common zero set must be zero according to equation 
©• So, we need to find those frequency components of 
the received codeword that are equal to those components 
of the error vector in the frequency domain. This results in 
solving fewer syndrome equations that finally gives us the 
error locations directly in the time domain. 

In ID codewords, the burst length is a code design parame¬ 
ter, and the shape of the burst is inconsequential. The situation 
is much more complicated in 2D. For example, a three bit burst 
error can be considered in three shapes as shown below: 



(c) 


Figure 1. (a) A burst error of three bits in horizontal shape, (b) A burst error 
of three bits in L-shape. (c) A burst error of three bits in vertical shape. 

The code is designed to correct specific error shapes as 
mentioned below. The same code construction is extended to 
correct unknown error patterns that are a disjoint combination 
of specific error patterns. We begin with a few definitions 
relating error patterns. 

Definition 1: A local 2D error pattern is defined as an 
error pattern which has a polynomial expression b(x, y) over 
GF( 2) such that the coefficient of the monomial x°y° in (x. y) 
is always 1. For a code of dimension n x m, the maximum 
degrees that x and y can take are n— 1 and m— 1 respectively. 

Definition 2: A global 2D error pattern is defined for a 
local error pattern b(x,y) starting from position (i,j). This is 
given by the following equation: 

e(x,y) = x l y J b{x,y). ( 20 ) 

As it is clear from Definition 1, an error pattern need not 
necessarily mean contiguous erroneous positions. This means 
that an isolated non-contiguous erroneous position is also 
considered as an error pattern. 

A. Predefined error patterns 

In this paper, we have considered horizontal and vertical 
error patterns as our basic error patterns. 


(a) 


X 

X X 


(b) 
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1) Error Pattern 1 (Horizontal burst): Let us consider a 
horizontal burst of error of dimension 1 x bi starting at position 
The error vector in the time domain is 

e\(x,y) = x l y^ + x l y^ +1 + • • • + x % y* +bl . ( 21 ) 


The 2D transformed error vector for (21 1 is given by 


n—1 m— 1 


E el = Y.Y. = 1*$* + 7 ie /3°' +1) ^ + 

i= 0 j= 0 

-h 'y ie /3 U+b (22) 





X 

X 




X 









(c) 


Figure 2. (a) Combination of a 1 x 2 and 2x1 error pattern forming a joint 
L-shaped burst, (b) Combination of a 1 x 2 and 2x1 error pattern forming 
two disjoint burst error, (c) Combination of a 1 x 2 and 2x1 error pattern 
overlapping at one position to form a smaller L-shaped burst. 


The following property holds true. 

Lemma 4: A horizontal burst of dimension 1 x (to — 1) i.e., 
putting L| = to, cannot be corrected. 

Proof: From \22\ , 


We restrict our code to correct errors with error pattern 1 or 
error pattern 2 or a combination of both. Let us assume these 
error patterns start at position and (£ 2 ,( 2 )- Consider 

the following error polynomial. 


43 = y e /3 i0 (i + /? 0 + /3 20 + --- + £ (bl_1)0 ), 

= 7 ie /3^(l + /? + /3 2 + ---+ ( 3 (bl - 1 >4 (23) 


e{x,y) = c 1 e 1 (x,y) + c 2 e 2 (x,y), 

= CiX kl y l1 (1 + j/) + c 2 x k2 y l2 (1 + at), (27) 


As f3 m = 1, we have, 

W + !)(! + P + !3 2 + /? 3 + ■ • • + /3 m_1 ) = 0. (24) 


where ci, C 2 £ GF( 2). ci = 1 and c 2 = 0 implies 
occurrence of error pattern 1 only, ci = 0 and C 2 = 1 implies 
occurrence of error pattern 2 only. Ci = 1 and C 2 = 1 implies 
occurrence of both error patterns. 


It is clear from the above equations that the maximum value 
that bi can take is to — 1. If b\ = to, from ( p4j >, 4'j, is always 
zero V( 0 , <f>) GVF C which indicates undetected error. ■ 
2) Error Pattern 2 (Vertical burst): Let us consider a 
vertical burst of error of dimension b 2 x 1 starting at position 
(i,j). The error vector in the time domain is 

e 2 (x, y) = x l yi + x l+1 y J + • • • + x z+b2 yF (25) 
The 2D transformed error vector is computed as 


7 i—l m— 1 

45 = 51 12 7*= 7 * V* + 7 (i+1) V i0 + 

2—0 j =0 

• • • + 7 ( i +k) fl J 0 J* (26) 

From Lemma 4 and using 7 " = 1, the maximum value b 2 
can take is n — 1 . 

From here on, we will consider b\ = 2 and b 2 = 2. 


B. Unknown error patterns 

Simultaneous occurrence of both horizontal and vertical 
burst errors could result in a number of unknown error 
patterns. Considering a horizontal error pattern of dimension 
1x2 and a vertical error pattern of dimension 2x1. The 
following figure shows some possible unknown error patterns. 


C. Selection of common zero set 

The common zero set specifies the code completely. The 
correction ability of the code is decided by selecting the 
common zero set. The common zero set is selected considering 
the following two salient points: 

• Disjoint syndrome sets among predefined error patterns. 

• Inclusion of indicator roots that help us to check the 
presence of a particular type of error. 

We elaborate these steps now. 

1) The disjoint burst criteria: According to the syn¬ 
drome components obtained at the common zero positions are 
equal to those corresponding components of the error vector 
in the frequency domain. We equate these set of equations and 
solve for i and j. A criteria for obtaining disjoint syndrome set 
is mentioned in If £-\ and E 2 are the set zeros of two error 
patterns bi(x,y) and b 2 (x,y) respectively, and if V c is the set 
of common zeros, then the following result must follow to get 
disjoint syndrome set between the aforesaid error patterns i.e., 

£1 n v c ± £2 n v c . (28) 


This can be understood in the frequency domain in a different 
way. Consider only horizontal and vertical error patterns. 
Equations ( |22] i and ( |26| ) give us a set of equations for solving 
the erroneous codeword for horizontal and vertical error pat¬ 
terns. When a horizontal burst error occurs, the set of equations 
given by ([22]) can be solved. On the other hand, the system 


of equations from (26 1 should result in no solution. If it does 
not happen, we have 

Jtepj<t> _|_ ye^(j+i)0 _ kopi<t> _|_ ^(k+i)e^14,^ ^9) 
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where, it is assumed that the syndrome for error of pattern 
1 occurring at position (i, j) matches with the syndrome for 
error of pattern 2 occurring at position (k 1 1). We have taken 
b \ = 2 and 62 = 2. Thus, for disjoint syndrome set, ( 29 1 
should have no solution among i, />;, and l. 

2) Indicator Zero set: 

Theorem 1 : Let Z b be the set of zeros of an error pattern 
b(x,y). Let V c be the set of common zeros for the code and 
V r be the roots of the received codeword. If the transmitted 
codeword c(x, y) is affected by b[x , y) and Z b r\V c 7 ^ 4> then 
V (A Xj A y ) £ V r 


(A X ,A y ) £ Zb n V c - (30) 


Proof: We have. 


This implies (X Xi X y ) £ Z b and (A Xj A y ) ^ V c . Thus, we 
have c(X X} X y ) = 0 which implies (A X) A y ) £ V c , which is 
again a contradiction. 

Case 3: (A Xi A y ) £ Z b n V c . 



r(x,y) = c(x, y) + b{x, y). 


(31) Figure 5. V r C (Z b nV c ). 


Let us consider (X x X y ) ^ Z b (T V c . 

This implies three following cases: 

Case 1, (X x X y 'j Z b and (A x A y ) ^ V c . 



Figure 3. V r D (V c U Z b ) = <f>. 

Setting (x,y) = (X x> X y ) we get, 

f (A Xj A y) = C (A X Xy) + b (X X ^Xy) . (32) 

Clearly the R.H.S is not equal to zero which shows (A Xi A y ) 
is not a root of r(x,y), which is a clear contradiction. 

Case 2: (A Xj A y ) £ Z b n V c . 



Figure 4. V r C (Z b n W) ■ 


This implies (A Xi A y ) £ V c and {X x X y ) ^ Z b . Thus, we 
have c(X x X y ) = 0 which implies {X x X y ) £ V c which is also 
a contradiction. 

Thus, we can conclude V (X x .X y ) £ V r , 

(A X ,A y) £ Z b D V c . 

■ 

Definition 3: Indicator zero set is a subset of the common 
zero set which is used to identify the presence of an error 
pattern a-priori , but not its locations. Each element of the 
indicator zero set is also a root of at least one of the error 
pattern. 

Theorem 1 shows the necessity of including at least one 
zero from each possible error pattern into the common zero 
set. These set of zeros form the indicator zero set. Evaluating 
the received codeword at these zero values help us to identify 
the error pattern. The other zeros in the common zero set 
are not chosen from any set of roots of the predefined error 
patterns. The syndromes at these positions will be non-zero if 
the received codeword is erroneous. The non-zero syndromes 
are used to construct the syndrome equation which will finally 
be used to solve the error locations. 

Theorem 2: Let Z b be the set of zeros for the error pattern 
b(x, y), and let V r be the set of zeros for the received codeword 
r{x. y). If Z b = then it cannot identified a-priori and V r = 
fi. 

Proof: If c(x, y ) is the transmitted codeword affected 
by b(x,y), the received codeword r{x,y) is given by the 
following equation: 

r(x,y) = c(x, y) +b(x, y). 

For all possible values ( x,y ), we have b(x,y) 7 ^ 0. Thus, 
r(x,y) 7 ^ 0, which means V, = </>. Thus, we cannot identify 
such error patterns, as there could be more than one error 
pattern having empty root sets. ■ 

In this paper the predefined error patterns considered are of 
dimension 1 x 2 and 2 x 1 . 
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The bi-variate polynomial expression for the 2D local error 
patterns are, 

bi(x,y) = 1 + y, (33) 

b 2 (x,y) = 1 + x. (34) 

Consider a code of dimension n x rn. The respective set of 


roots for the above 

equations in time domain are. 


4 1} = 

{(1)1)! (7>1)) ' ' 


(35) 

7 (2) 



(36) 


Presence of both error patterns starting from different lo¬ 
cations and not overlapping has the following polynomial 
equation. 


e 3 (x,y) = x kl y h (l + y)+x k 2 y l2 (l + x). (37) 

The roots of this equation are given by 

zf ={(1,1)}. (38) 

Thus, from Theorem 1, the common zero set should have 
the root (1,1) to identify the presence of e 3 (x, y). To identify 
ei (x,y) and e 2 (x,y), we have to select a root from Z^ and 
Z^ respectively. The following set vi^ gives us the set of 
indicator zeros contained in V c . 

V< J) ={( 1 , 1 ), ( 7 ,1) ,(!,/?)}■ (39) 

Thus, when r(x, y) is received, we evaluate it at the roots from 
the set Vc'' and identify the presence of an error pattern. We 
have the following cases: 

Case 1: e^ix^y) is identified as the error pattern if, 

. r(l, 1 ) = 0 . 

• r( 7 , 1 ) ^ 0 . 

. r(l,/3)^0. 

Case 2: e\(x,y) is identified as the error pattern if, 

. r(l, 1 ) = 0 . 

• r ( 7 , 1 ) = 0 . 

. r( 1 , 0 ) ^ 0 . 

Case 3: e 2 (x,y) is identified as the error pattern if, 

. r(l, 1 ) = 0 . 

• r( 7 , 1 ) ^ 0 . 

• r (1, 0) = 0. 

Case 4: No error pattern is identified if, 

. r(l, 1 )^ 0 , r (7, 1)^0 and r ( 1 , 0) ^ 0 . 

• r(l, 1 ) = 0 , r (7, 1 ) = 0 and r ( 1 , 0) = 0 . 

• r(l, 1 )^ 0 , r (7, 1 ) =/= 0 and r ( 1 , 0) = 0 . 

• r(l, 1 ) ^ 0 , r (7, 1 ) = 0 and r ( 1 , 0) ^ 0 . 

• r(l, 1 )^ 0 , r (7, 1 ) = 0 and r ( 1 , 0) = 0 . 

The last case signifies a uncorrectable error pattern. We 
illustrate this through an example and formalize the procedure 
within Algorithm 1. 


Algorithm 1 Decoding Algorithm 
Input: received codeword r(x,y) and V c . 

Step 1: Indicator zeros. 

. Set vi J) ■<— {(1,1), ( 7 , 1), (1, 0)} ■ 

. if r(l, 1) = 0 && r( 7 ,1) ^ 0 && r(l, /3) ^ 0 

Set ci 1. 

Set c 2 <— 1. 

. else if r(l, 1) = 0 && r( 7 ,1) = 0 && r(l, /3) ^ 0 
Set ci 1. 

Set c 2 <r- 0. 

. else if r(l, 1) = 0 && r( 7 , 1) ^ 0 && r(l, /3) = 0 
Set ci 0. 

Set c 2 <r- 1. 

. else if r(l, 1) = 0 && r( 7 , 1) = 0 && r(l, /3) = 0 
Proceed to Step 2. 

• else 

Output “Error cannot be corrected”. 

Step 2: Solve syndrome equations 

• Evaluate syndrome values using the remaining values in 
V c on r(x, y). 

• if all syndrome values are zero. 

Set c(x,y) <- r(x,y). 

Go to Output. 

• else 

Solve the syndrome equations. 

Flip the solved erroneous positions to get c(x,y). 

Go to Output. 

Output: c(x,y). 


Example 2: Let us assume that we have transmitted the zero 
codeword. Our basic error patterns are 1 x 2 and 2x1. Let 
the common zero set in the frequency domain be, 

V c = {(1,1)( 7 ,1),(1,/3),(7,/?),(7 2 ,/3 3 )}- (40) 

Consider the following two received codewords. 


/ 1 

1 

1 

0 

0 \ 

( 0 

0 

0 

0 

0 \ 

-i 

II 

0 

0 

1 

0 

0 

r 2 = 0 

0 

1 

1 

0 

V 0 

0 

0 

0 

0 ) 

V 0 

0 

0 

0 

0 ) 


The syndrome equation considered here is given by the 
following equation. 


EgJ = Cl ^7* 1 <y i * + 7 * 1 *0 ( « 1 + 1) *) + 

c 2 +7 (fc 2 +i _ ( 41 ) 

The elements (0,0), (0,1) and (1, 0) are the indicator zeros. 
Decoding 77 : We have, 

n(x,y) = 1 +y + y 2 +xy 2 . (42) 

It is clear that 77 ( 1 ,1) = 0, ri( 7 ,1) = a 10 and 
ri(l,/3) = a 14 . Thus, the error polynomial is of the 



form equivalent to (x,y). Hence, we make C\ = 1 and 
C2 = 1. To solve for the error locations, we plug in 

(7,1), (1, ft), (7, ft), and (7 2 , ft 3 ) in equation (371 to obtain 
the following syndrome equations: 


Theorem 3: Presence of 7 number of burst errors of type 
1 x A can be uniquely identified. 

Proof: Presence of // number of burst errors of type 1 x A 
leads to the following expression in the transformed domain: 


v fc 2 


. ft( l i+l) 


II 

0 

M-l /A—1 

5Z (H i ia0 ft [ ' ja 

14 

— O' 

0 

II 

3 

0 

II 

b 

= « 7 ,’ 

Equation ([46) can be simplified 


M-1 

= a . 

M 

r 

q 

II 

(0,0) and 

(7 = 0 


as 






(46) 


(47) 


(k 2 l 2 ) = (0,2). 

Decoding r 2 : We have, 

n{x,y) = xy 2 +xy 3 . (43) 

It is clear that ri(l,l) = 0, 77(7,1) = 0 and r\{l,ft) = 
a 5 . Thus, the error polynomial is of the form equivalent to 
ei(x,y). Hence, we make <7 = 1 and C 2 = 0. To solve for 
the error locations, we plug in (l,/3), (7, /3), and ("f 2 ,ft 3 ) in 
ei (x,y) = x kl y ll ( 1 + y) to obtain the following syndrome 
equations: 

fth + ft(h+l) = a 5, 

7 fc ijS‘i+7 fc i/9 (,1+ i) = a 10 , 

7 2 fci^3ii _|_ ^,2fcie^3(ii+l) _ Q 5 


Let {(Ccr, »7cr)} and {{Ca,^)} for 0 < a < p - 1 be two 
solutions of the above equation with none of component of the 
solutions equal. As the right hand side of the above equation 
does not have a variable, we have 

fl— 1 (JL — 1 

J^^^ft^ + ^2^/3^ = 0. (48) 

a —0 <7=0 

Equation ( |48] > implies either the terms are all equal, or they 
are pairwise equal. This gives us the following ordered pairs: 


(Co, ^o) = 

(C o, Vo), 


(Ci, 7i) = 

(Ci , Vi), 


(Cm,^) = 

(c;,<). 

(49) 


Solving the above equation, we have (k\l\) = (1,2). 
Proposition!: Presence of overlapping predefined error pat¬ 
terns of dimension 1x2 and 2x1 cannot be determined 
a-priori and hence, cannot be corrected. 

Proof: The polynomial expression of the error pattern as 
a linear combination of the error pattern of dimension 1x2 
and 1 x 2 at position (fci.Zi) and (^2,(2) respectively is given 

by, 

e(x, y) = x kl y l1 (1 + y) + x k2 y l2 (1 + x). (44) 

Without loss of generality, let k 2 = k\ and l 2 = l\ + 1 for 
overlapping errors. We get, 

e(x, y) = x kl y h + x kl+1 y h+1 , (45) 


Any other combination will just be a permutation of the 
solution set which means, the solution set is unique, proving 
the result. ■ 

The following example will illustrate our idea. 

Example 3: Let the received codeword be 

0 0 0 1 1 \ 

0 0 0 0 0. 

1 1 0 0 0 / 

Let the two bursts start at positions (i, j) and (k,l). The 
corresponding error polynomial is 

E e j = ^ftrt + 7 *ftU+V* + ^ ft 1 * + (50) 

Let the common zero in the frequency domain be 


which gives two isolated errors. Such an overlapping burst can 
be treated as a disjoint combination of an error of dimension 
1x2 error and an isolated error, or as a disjoint combination 
of an error of dimension 2x1 and an isolated error. From 
Theorem 2, we have established that an isolated error is not 
within our set of predefined error patterns. Thus, an error 
pattern formed in combination of any valid predefined error 
pattern with an isolated error cannot be identified a-priori. ■ 

D. Multiple Bursts of Same Type 


= {(0,0), (1,1), (1,4), (2,2), (2,3)}. (51) 

This gives us the following equations as before i.e., 

-E' 0 ,o = 0 = 7?o,o, 

E ltl = ^fti + 7 i /3(7+ 1 ) + ^kpi + 7 *0(i+1) = a i2 ) 

E 1A = Yft^ + 7 i /3 4( 4+ 1 > + y k ft 41 + 7 fc /3 4 ('+ 1 ) = a 3 , 

E 2 ,2 = l 2l ft 2j + 7 2i /3 2(j+1) + y 2k ft 21 + 7 2fc / 3 2 ( ; + 1 ) = a 9 , 
E 2 , 3 = l 2i ft 3j + 7 2i ftHj+l) + rfkpZl + 7 2fc^3(i+l) = a 6_ 


If the total number of bits in error is within the minimum 
distance bound, this code can be used to correct the presence 
of multiple bursts of an error of a particular type. Let us 
assume that the bursts are non-overlapping. For the general 
case, suppose there are p number of burst errors of type 1 x A 
that are non-overlapping, we can uniquely identify the position 
of the errors. The following theorem summarizes this result. 


The solution to the above system of equations gives the 
result (i,j) = (0,3) and (k,l) = (2,0) which are exactly the 
positions where the 1x2 bursts start. 

Previous decoding schemes fH required storing the syn¬ 
drome polynomials for each error pattern. In our case, the 
erroneous positions are found directly by solving the appro¬ 
priate syndrome equations. 
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V. Conclusions 

We presented a 2D cyclic code construction in the frequency 
domain. Our code is capable of correcting simultaneous oc¬ 
currence of multiple predefined error patterns. The encoding 
and decoding steps are done entirely in the transformed 
domain, thereby, facilitating efficient realization. Further work 
is needed to extend this to a broad framework of 2D codes 
over curves. 
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